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Abstract 

We study momentum relaxation due to dilute, weak impurities in a strongly coupled CFT, 
a truncation of the M2 brane theory. Using the AdS/CFT correspondence, we compute the 
relaxation timescale as a function of the background magnetic field B and charge density p. 
The theory admits two different types of impurities. We find that for magnetic impurities, 
momentum relaxation due to the impurity is suppressed by a background B or p. For electric 
impurities, due to an underlying instability in the theory towards an ordered phase, the in- 
verse relaxation timescale increases dramatically near ^B 2 + p 2 /a 2 ~ 21T 2 . We compute the 
Nernst response for the impure theory, and comment on similarities with recent measurements 
in superconductors. 



1 Introduction 



Second order quantum phase transitions can occur at zero temperature in condensed matter systems 
when there is a non-analytic change in the ground state energy as a function of some coupling [1 . On 
either side of the critical point, at sufficiently low temperatures the system admits a quasiparticlc 
description. The nature of the quasiparticles depends on the orders characterising the different 
phases. For instance, Neel ordered phases have spin wave excitations whereas phases with Valence 
Bond Solid order have spinon and 'photon' excitations [2j. However, as the temperature is increased, 
near to the critical coupling, thermal fluctuations lead to competition between the different orders. 
In this strongly coupled regime, in which both quantum and thermal fluctuations play an important 
role, neither quasiparticle description is appropriate. Instead the system is best described by a finite 
temperature 2+1 dimensional conformal field theory (CFT) [JJ[2]. Often this theory turns out to be 
relativistically invariant (the 'speed of light' in these theories is not c but rather some lower speed 
characterising the material). 

Quantum critical points are believed to be important in various systems of experimental interest, 
including the high T c cuprate superconductors [3J3]. In contrast, the theorist's toolkit of tractable 
2+1 CFTs has largely been limited to the O(N) model at large N, and related models with relevant 
quartic interactions [TJ [2] . The AdS / CFT correspondence [5] provides a wealth of new examples of 
2+1 CFTs in which explicit computations are possible at large N. Although this correspondence 
has been intensely studied for a decade, potential applications to concrete condensed matter systems 
have only recently begun to be explored [5J [7J El [5] . Relevant earlier work on finite temperature 
physics in AdS/± includes [lOj [Til H2 E] . Most immediately these CFTs provide new solvable toy 
models for strongly coupled dynamics in 2+1 dimensions. In the future it might be possible to 
engineer the relevant supersymmetric CFTs in a lab. One can thus speculate that in addition to 
having interesting physical properties in their own right, these systems open up the possibility of 
doing experimental quantum gravity via the AdS/CFT correspondence. 

In many experimental systems, microscopic impurities in the samples leave a significant imprint 
on the physics. Most directly, the impurities are at fixed (random) locations and therefore break 
translational invariance. This loss leads to a relaxation of momentum at late times. Momentum 
relaxation is characterised by a timescale r; mp , which will be the main focus of this paper. In many 
contexts, see e.g. [8] for a discussion, momentum relaxation is necessary for dc (oj — 0) transport 
to be finite; otherwise the 'Drude peak' becomes a delta function as 1/Tj mp — > 0. We will exhibit 
similar 'metallic' behaviour in a strongly coupled CFT with impurities. 

Given a theory, the effect of impurities can in principle be computed microscopically. One would 
like to know how the impurity timescale depends on temperature, background magnetic field and 
charge density. To our knowledge there is no first principles computation of 7i mp (T, B, p) available 
in the literature. If the impurities are sufficiently dilute and furthermore weak enough for their 
interaction with the CFT to be treated perturbatively, then rj mp can be computed from a two point 
function in the finite temperature CFT. We will derive the formula for T; mp and use the AdS/CFT 
correspondence to compute it for a truncation of the M2 brane theory. 

In sections 2 and 3 we give a general description of the effect of dilute, weak impurities in a CFT. 
The main result of these sections is equation ( 12 1, which expresses the impurity relaxation timescale 
in terms of the retarded Green's function of a relevant scalar operator in the theory which couples 
to the impurities. 

In section 4, we show how this Green's function may be evaluated in a truncation of the M2 
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brane theory using the AdS/CFT correspondence. Our truncation is not consistent. We neglect 
a coupling between the scalar and pseudoscalar modes we keep and a second gauge field. The 
inconsistent truncation is made to keep the calculations simple and also because we have no good 
physical interpretation of the second gauge field in our condensed matter context. In another small 
abuse of the M2 brane theory, we give our pseudoscalar operator conformal dimension one instead of 
two. This gives us a larger collection of relevant operators to study. Note however that the standard 
M2 brane theory is obtainable as a double trace deformation of ours [H] . 

In the following section 5 we give our (numerical) results for T; mp (T, B, p). We first note that 
there are in fact two relevant operators in the theory, dual to linear combinations of the scalar and 
pseudoscalar of the bulk, with considerably different effects. The operators can be distinguished by 
their transformation properties under charge conjugation C, parity V, and time reversal T. In the 
presence of a background magnetic field and finite charge density, the first operator, C + , transforms 
the same way as the magnetic field while 0_ transforms as the electric field. We are thus tempted 
to identify 0+ with magnetic impurities and 0_ with electric impurities. 

For both operators, we find that the impurity time l/ri mp depends on the charge density p and 
magnetic field B only in the combination \J B 2 + p 2 Joq where ctq is the electrical conductivity at 
p = B = 0. This combined dependence, while remarkable, may be an accident of the details of the 
M2 brane theory. 

For the 0+ impurities, we find that increasing the magnetic field or charge density suppresses 
momentum relaxation due to impurities. At the same time, at nonzero B there is an independent 
source of momentum relaxation due to hydrodynamic cyclotron motion. Increasing B increases the 
momentum relaxation due to cyclotron motion. 

However, for the other operator, 0_, momentum relaxation becomes faster with increasing B or 
p. When \J B 2 + p 2 /erg reaches a critical value of about 21T 2 , then l/ri mp — > oo. We will see that 
this divergence is symptomatic of an underlying instability of the theory towards an ordered phase. 
From a gravitational point of view, the instability occurs when the scalar field dual to 0_ develops 
an unstable mode. This instability provides a clean counterexample to the original version of the 
Gubser-Mitra 'correlated stability conjecture' [TBI [T7] , similar to those discussed in [TFj . 

In section 6 we use our results to compute the Nernst coefficient in the CFT. We note that our 
plots show some qualitative similarities to the organic superconductors studied in [19], perhaps indi- 
cating a nearby quantum critical point in that system. While the organic superconductors effectively 
have vanishing charge density, p = 0, because they are not doped, in the cuprate superconductors 
the doping x can be thought of as a nonzero charge density [5] . 

2 The effect of random disorder 

The presence of impurities breaks the translational and conformal invariance of the CFT. Locally 
the impurities will source relevant operators and potentially drive the theory away from the fixed 
point. Thus we can model the impurities by adding the following coupling to the Hamiltonian 



The operator 0(t,y) is the most relevant operator in the conformal field theory that preserves the 
global symmetries of the theory. Charged impurities are also of interest, although we shall only 
study neutral impurities here. This term breaks translation invariance, because the potential V(y) 




(1) 
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is explicitly space dependent, and is hard to work with directly. To get around this one treats 
the impurity potential statistically. If the impurities are sufficiently dilute that their effects do not 
overlap, then the precise weighting on the space of potentials is not important, and it is useful to 
take the Gaussian J20J 

(...) imp = fcVf/W^f...), (2) 



which implies 

<n*)>im P = , (V(x)V(y)) imp = V 2 5 ( - 2 \x - y) . (3) 

We will furthermore make the approximation that scattering off impurities may be treated pertur- 
batively, so that we can expand in powers of V . The strength of the potential V is a dimensionful 
quantity, with mass dimension 

[V]=2-A a . (4) 

So long as 2 — Aq > 0, V has positive scaling dimension and the impurities arc a relevant perturba- 
tion. This condition that V has positive scaling dimension is often called the Harris criterion [2T] , 
In the case of the M2 brane theory, the most relevant neutral operators are mass terms for the scalar 
fields [22] [23j [24] . By 'neutral' we mean with respect to a certain U(l) C 50(8), as we describe 
below. These operators have conformal dimension Aq = 1. We can anticipate that the relevance of 
these operators implies that they will become strongly coupled at low energies. We will be careful 
to keep track of the regime in which perturbation theory in V is valid. 

The breaking of translation invariance leads to the late time non conservation of momentum. 
The timescale associated with the loss of momentum is called the impurity relaxation timescale and 
is given to leading order in strength of the impurity potential, V, by 

J_ = _J, lim ^M^). (5) 

Timp X0 W 

This equation follows from considering the 'memory function'. We give a derivation in the following 
section, along with a discussion of the precise meaning of Ti mp and of regimes of validity. To fix 
conventions, we define the retarded Green's function of any operator O in Fourier space to be 

oo 



G§ (w,k) = -i J d 2 x J dt([O(t,x),O(0,0)})e wt - lk - x . (6) 

The quantities that appear in §5§ are firstly 

Xo =KmGR p (Lj,0) = e + P, (7) 

where V is the momentum density in a fixed direction, i.e. V = riiT 01 for some unit vector n £ M 2 . 
The number xo is the static susceptibility for the momentum density, e is the energy density and P 
the pressure. The remaining term in ^ is the Green's function for the operator 

T{t,x) = [P(t,x),SH] = J d 2 yV(y)[V(t,x),0(t,y)] . (8) 

This expression simplifies using that the spatial integral of V(t, x) is a momentum and hence gener- 
ates translations. Locally therefore [V(t,x),0(t,y)] — i8^ 2 \x — y)dO{t,y) up to a total derivative 
with respect to x. We shall ignore this total derivative; one can use translation invariance to show 
that this derivative does not contribute to the Green's function Gjr-p. Thus ^ becomes 

J 7 (t,x)=i / d 2 yV(y)S (2) (x - y)dO{t,y) = iV(x)dO(t,x) . (9) 
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Here d denotes spatial derivative in the same direction as V, i.e. d = n % di 
The Green's function is therefore 



G%r(u,0) = iV 2 dt([dO(t,0),dO(0,0)])e lLUt , (10) 



where we used ([3]). Passing to momentum space, one obtains 

?2 r d 2 k 
2 J (2ny 

The factor of 1/2 arises because the integral with (n • fc) 2 in the integrand is half the integral with 
fc 2 , by isotropy. Thus 



G^( W| 0) = --V / j— fc 2 Gg (c,fc). (11) 



1 - ^ l im /• d 2 k fe2 ImGg (a;,fc) 



' imp 



2x0^^0 7 (27r) 2 u; 



(12) 



This is the formula we will use to compute the relaxation timescale T; mp . We will obtain the 
dependence on charge density p and background magnetic field B. Due to conformal invariance and 
dimensional analysis, one has the scaling form 

1 V 2 ( p B 



r imp t*-**° f \t 2 't 2 ) ■ (13) 

3 The impurity timescale 

3.1 Impurity in the absence of a magnetic field 

In this section we derive the expression ([5| for the impurity timescale. First note that from the 
definition of the retarded Green's function, and (twice) using the Heisenberg equation of motion one 
obtains 

uj 2 G* v (w, fc) = -G$ >HWtH] (uj, k) + Gf V H][V H] (0, k) . (14) 

This equation is true up to contact terms for any operator. In the absence of the impurity potential, V 
is the conserved momentum density as well as the energy current. If we split the total Hamiltonian 
as H = Hq + SH, with SH given by ([lj above, then Hq is the Hamiltonian of a translationally 
invariant theory while SH explicitly breaks this symmetry. It follows that [j d 2 x V(x), Hq] = 0. 
Thus from the definition of the retarded Green's function, at zero spatial momentum we have that 
G^, Ho ]o(u, 0) = 0, for any operator O. It follows that 

^ 2 G^,0) = -(G^( W ,0)-G^(0,0)) . (15) 

Evaluated in a background with a nonzero V, this result for G^-p is exact. For us the usefulness 
of this expression lies in the fact that a V 2 trivially factors out of the right hand side, see equation 



( 10 ), allowing us to evaluate the remainder in a background with vanishing V , thus yielding a result 



for Gp.p that is accurate to leading order in V 2 . 



Given this formula (151, limy^ G^p V (ui 1 0) = as required by the restoration of transnational in- 
variance in this limit. However, we have no guarantee that at fixed nonzero V the lim^^o G^-p(uj, 0)/u>, 
which is proportional to a dc thermal conductivity, is finite. Yet physically, we expect the system 
with impurities to be well behaved in response to low frequency perturbations. These considerations 
suggest that the response of the system is governed by a modified Green's function where a contact 
term has been subtracted: 

i9(w,0)Tw = G^(u,0)-xo. (16) 
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The quantity 9(uj,0) might be called a momentum conductivity; in appendix [A] we show that indeed 
(V) = — #(cj,0)VT. By 9 we mean a diagonal component of the conductivity matrix; there is no 
off-diagonal conduction in the absence of a magnetic field. 



We can now read off xo m ( 16 )■ In a translationally invariant system in the absence of a magnetic 
field, the momentum conductivity #(w,0) = i(e + P)/loT where e is the energy density and P the 
pressure. This result can be seen from hydrodynamics [5] or from the Ward identity arguments in 
[9] . The fact that the imaginary part of the conductivity behaves like 1 jio implies that the real part 
will contain a 5(u>), because with some suitable ie prescription l/oj — P(\/lo) — in8(uj). This delta 
function is expected for the energy flow in a translationally invariant system, which has no way to 
dissipate momentum. From the fact that translation invariance implies Gpp(w,0) = 0, it follows 
that 

Xo = ( + P. (17) 

That xo = lini^^o G^-p(uj, 0) remains finite as V — ► is consistent with the fact that the lj — > and 
V — > limits of G pv (cj, 0) do not commute^] 

The memory function method 

Now we extract the timescale from G^ v (lu, 0). The memory function M(ui) is defined by 



uG* p (u,0) 



which can be formally inverted to yield 



G -^°) = ^TT- < 19 > 

LO + M (LO) 

The reason for introducing this function is to reinterpret the small lo limit of M(w) as an inverse 
scattering time i/Ti mp . We comment on this interpretation below. At finite lo, by translational 
invariance we expect G pv {lo,Q) to vanish in the absence of impurities. We also expect that the 



inverse scattering time goes to zero with no impurities. Thus, consistent with eq. (18 1, we take both 
M(lo) and G pv (uj,0) to scale as V 2 . These statements do not hold when there is a nonvanishing 
background magnetic field. For the moment we are setting B — 0. As we noted, xo = G^-p (0,0) 



is independent of V to leading order. Thus from eq. (18 1, and using equation (15 1, we find the 
approximate expression for M(lo): 

^ uG$ v {u,0) = -(G^,0)-G^(0,0))/a; ^ 
Xo Xo 

The approximation is up to terms of 0(V 3 ). We would like to take the small to limit of this expression. 
Because V is dimensionful, one should expect that in order to obtain the true to — > limit it will be 
necessary to resum higher order contributions in V /ui 2 ~ A ° . Indeed, we have already noted that the 



lo — ► and limits do not commute. Using (20) directly to compute the relaxation timescale, 



is often called the 'memory function method', see for instance refs. [25, 2B]. It was shown in ref. [57] 
that this method generally does not give the correct answer for small frequencies lo 2 ~ a ° < V. 

We can obtain reliable results in the case when the temperature is large compared to the strength 
of the scattering potential: V" 1 ^ 2_A °- ) -C T. Given this separation of scales, at low frequencies 



1 In the absence of an external electric field, 9 is related to the thermal conductivity R and thermoelectric coefficient 
a defined, for example, in ref. [9] via 8 = R + &fi. Thus at vanishing chemical potential fi, the thermal conductivity 
is 8. With no chemical potential, then xo = sT, with s the entropy density, leading to R = is/w. 
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to <C T, the constraint lo 2 ~ a ° > V, for the validity of perturbation theory, should be replaced by 
the weaker constraint T 2 ~ A ° > V. In the language of the renormalization group, V is a relevant 
operator and we must choose a scale at which to evaluate it. At nonzero temperature and at 
frequencies lo < T, we expect the temperature to act as a cut off in the renormalization flow for 
V. In higher order corrections to the memory function M, the ratio V/T 2 - A ° will appear in place 
of V /lo 2 ~ a ° . Thus at high temperatures, we are justified in treating the impurities perturbatively, 
including at very low frequencies. 



Evaluated perturbatively, the expression for the memory function ( 20 ) has an overall factor of 



V 2 , and no other dependence on V. Scaling therefore implies that the frequency dependence must 
be some function of lo/T. Thus the hydrodynamic limit lo ^ T is equivalent to taking lo — * in 

_(G^( w ,0)-G^(0,0)) 



( 20 ) . We can identify 



' imp 



— lim 

Xo o 

1 dG^(u,0) 



(21) 



Xo 



duo 



Using the fact that Re G^-p is an even function while Im G^-p is an odd function of lo, we 
main result of the section that 



find 



1 ,. ImG^(w,0) 



= lim 



(22) 



Let us clarify the physical meaning of Tj mp . Strictly, the momentum relaxation timescale is given 
by the imaginary part of the pole in G^ v (lo, 0) that is closest to the real axis. The identification of 
the zero frequency limit as the inverse relaxation timescale 



lim M(lo) 



(23) 



is valid so long as the putative pole at — i/Vi mp is sufficiently close to the real axis, in the sense that 

1 



M 



(24) 



It is easy to check, for instance, that this will be true if G^ v (lo, 0) has several poles, but one is much 
closer to the real axis than the others. Given that 1/Tj mp ~ V 2 /T 3 ~ 2A ° <C T is parametrically small 



compared to the temperature scale, we expect that (24 1 holds. However, if the coefficient of V in 



( 12 ) becomes sufficiently large, then that expression is no longer reliable as a relaxation timescale. 

We argued above that at high temperatures, potentially dangerous V/lo 2 ~ a ° corrections would 
be replaced by corrections in V jT 2 ~ A ° . While V/lo 2 ~ a ° corrections would tend to add more poles 
close to the origin of the complex lo plane, V/T 2 ~ A ° corrections are suppressed. 



3.2 Impurity plus an external magnetic field 

There is great formal similarity between adding to the Lagrangian a random potential coupled to a 
scalar operator J d 3 x V(x)0(x) and adding an external U(l) gauge field coupled to a global current 
— J d 3 x A ll (x)J fl (x). This second case was considered in detail in ref. [3]. A constant magnetic 
field has F=dA=B dx A dy. Any potential A for this field strength will formally break translation 
invariance, which is one way to see why G^ v (lo,0) can be nonzero in the absence of impurities but 
in the presence of a background magnetic field. 



G 



From the Ward identity results in ref. [9], with V — one has 

uG^p(u, 0) = B 2 (G§j(cj, 0) - G%j(0, 0)) /u , (25) 

where J = Prii is the electric current. The uj —> limit of the current-current Green's function 
needs to be taken with care, as the B — ► and to — > limits do not commute. 

At V = 0, Gjjr(0,0) vanishes for real n*, but there is still a Hall conductivity [7J|S], 

lim gi^l^ = lp/B . (26) 

To see the Hall effect within our formalism, it is convenient to let the rii be complex and think of 
G R as an inner product. This method of complexifying the conductivities was found to be useful in 
[9]. In this paper we will do so for this section and appendix |b| only. Taking Uj = (1, — i)/y/2, we 
find 

— GjLjx — iG% Jy + iG JVJX + G% Jy . (27) 
By rotational invariance, this expression simplifies to 

Gjj — Cjxjx — iG% Jy . (28) 

Combining eqs. ( p5| and p^j ) then to include both impurities and a magnetic field leads to 

« 2 G£p(w, 0) = - (G$Au, 0) - G^(0, 0)) + B 2 {G*j(u;, 0) - up/B) . (29) 



Like eq. (151, we believe this expression is exact when the right hand side is evaluated in a background 



with nontrivial B and V and Gj = 0. We have been loose in our derivation, but the result can be 
made rigorous through a Ward identity argument. The potential cross term proportional to Gj 
will vanish anyway to leading order in V after averaging over V. But in fact, we will see shortly 
that the Green's function Gj = vanishes for the M2 brane theory at leading order in l/N before 
averaging. 

We call poles hydrodynamic if they are close to the origin of the complex u> plane compared to 
the scales T and p. In the previous section, in the absence of a B field, we argued that for small 
V there was a hydrodynamic pole in G-p-p at — ?/ri mp . In ref. [5], we saw that in the presence of a 
small magnetic field but in the absence of impurities, there are a pair of hydrodynamic poles at 

. B (±p - ioB) 

±u c -vy= — , 30 

e + P 

both in G-p-p and Gjj. These poles correspond to damped relativistic cyclotron motiorj^J Here a 
is the electrical conductivity of the CFT. 

Physically, if we have a small but nonzero V or B, we expect the u> — > response of the system 
to be well behaved. The poles in the Green's functions should interpolate continuously between the 



cases in which cither V or B vanish. Thus, contemplating (29 1, it is natural to expect that the 
hydrodynamic poles in the presence of both V and B are 

~ "iip - *7 ± w c . (31) 
This combination of T- lmp and the cyclotron pole was also found in the hydrodynamic analysis of [5] . 



In eq. (31) we see that there are two sources for momentum relaxation; one is impurities and the 



2 Again, our time dependence is e luJt . Thus 7 > and Ti mp > lead to exponential damping at long times. 
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other can naively be thought of as due to collisions of positively and negatively charged excitations 
of the theory undergoing cyclotron motion in opposite directions. In appendix [5] we demonstrate 
how the memory function method can be used to extract the correct poles from the hydrodynamic 
expressions for G-p-p and Gjj, and we also show that these hydrodynamic expressions satisfy the 
relation ( 29 1 . 



The impurity relaxation timescale Ti mp appearing in eq. (31) is given by the same formula as 
before, eq. (22 1, but can now have a B dependence, provided B is kept small so that our expression 
for the cyclotron pole is reliable. In more detail, our hydrodynamic expression for 7 is correct to 
order B 2 but we expect it to receive unknown corrections at 0(V 2 ). This is because the e + P 
appearing in the definition of 7 in (30 1 will have a V 2 dependence. Our expression (22 1 for 1/r 
is correct to order V 2 and in principle we know its full B dependence. However, we expect the 
hydrodynamic approximation to fail at 0(B 3 ). Thus, 1/r should be trustworthy at order 0(B 2 V 2 ) 
but no higher. Problematically, we do not know the corresponding 0(B 2 V 2 ) correction to 7. This 
uncertainty turns out not to be important for many of our considerations. In particular, when 
studying the Nernst effect at p = we will find that there is no dependence on 7. 



4 Impurity relaxation in the truncated M2 brane theory 

We have shown that computing the impurity relaxation timescale boils down to the two point 



function of equation (12 1. In this section we will perform the calculation for a truncation of the 
M2 brane CFT, using an adaptation of the AdS/CFT dictionary for real time two point functions 
proposed in [55]. As noted in the introduction, the truncation of the M2 brane theory is not 
consistent, and we make a nonstandard choice of conformal dimension for one of the scalars. 

The CFT at finite temperature with a nonzero charge density and background magnetic field is 
dual to a dyonic black hole in AdSi |7J. Such a black hole is a solution to Einstein-Maxwell theory 
with a negative cosmological constant. To this gravitational theory we would like to add a neutral 
scalar field (j) dual to our relevant operator O. 

From the M2 brane point of view, this Einstein-Maxwell theory with a neutral scalar should be 
a sector of eleven dimensional supergravity compactified on an S 7 . After the compactification, we 
are left with an Af = 8 supersymmetric gravitational theory on a four dimensional spacetime with 
a negative cosmological constant and infinite towers of Kaluza-Klein states. The isometry group 
SO (8) of the S 7 determines many features of the lower dimensional theory. 

It is believed to be consistent to truncate this compactification to the lowest states in the Kaluza- 
Klein towers. This truncation is called Af — 8 gauged supergravity. The bosonic content of the 
truncation is the four dimensional graviton, an SO(8) gauge field, and a scalar and pseudoscalar 
transforming under 35 dimensional representations of 5*0(8). These fields and their superpartners 
are thought not to act as sources for the other excited modes in the Kaluza-Klein towers, and hence 
the truncation is called consistent. 

There exists a further truncation of Af = 8 gauged supergravity to an abelian sector [29 . Consider 
the Cartan subalgebra of so(8). The group SO(8) acts naturally on M 8 , and we can think of the 
elements of the Cartan subalgebra as generators of rotations in the 12, 34, 56, and 78 planes of 
M 8 . The action for the abelian truncation involves only the four U(l) gauge fields corresponding to 
these four rotations and the three scalar and pseudoscalar fields neutral under these four U(l) gauge 
groups. The full action for this abelian truncation can be found for example in ref. |30j . 

To see which scalar fields remain neutral, it is useful to think of the 35 dimensional representation 
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of 5*0(8) as the set of symmetric traceless quadratic polynomials in the coordinates on R 8 , x % xK 
The ?7(l)'s generated by the Cartan subalgebra act as phase rotations on the complex combinations 
Zi = x 2 " 1 ^ 1 ±ix 21 for i = 1,2,3,4. Thus the neutral quadratic scalars must be \zi\ 2 — |zi + i| 2 . There 
are three independent such scalars. By allowing the dual operators O to couple to the impurity 
potential, we are breaking the SO(8) symmetry of the M2 brane theory. However, the operators 
are neutral under the U(l) symmetry for which we have a background B field and charge density. 
Indeed, models for 'real world' electronic systems often involve neutral impurities (see for example 
[ST]). 

We will make a further truncation of this theory and consider only the gauge field corresponding to 
a simultaneous rotation in all four planes, along with the scalar and pseudoscalar pair corresponding 
to just one of these neutral elements of the 35 dimensional representation. As we consider scalar 
two-point functions, we will only need the action to quadratic order in the scalars: 

\ f d^x^-g (R - | {(d^)(d^) + (d^xWx)} + 2£- 2 [3 + \(<t> 2 + x 2 )] 

-L 2 [1 + |(0 2 - x 2 )] F IIV F ,XV + iL 2 4>xe^F^F p(T ) . (32) 



S = - 

2k 



Here d tvpa is the totally antisymmetric tensor with e 0123 = 1/ y/^g, and L is a length scale which 
determines the AdS radius. The dyonic black hole is a solution to the equations of motion provided 



that 4> — X — 0- As we can see in (32 1, however, the fact that e tJ,vpa F^ v F pa does not vanish for the 
dyonic background implies that fluctuations in <fi and x source each other. This mixing is why we 
keep the two scalar fields. We comment below on the physical implications of the doubling. 

We note that in a charged background, our final truncation to the two scalars plus the diagonal 
Maxwell field is not consistent. We have neglected a quadratic coupling between the neutral scalar 
fields and fluctuations of a non-diagonal combination of the four U(l) gauge fields of the full M2 
brane theory [30] . Specifically, we are throwing away the terms 

S G = - A / d 4 x^g~ {—I?G liV G tLV - 2L 2 <j>F^G^ + L\e^° 'F^G pa ) , (33) 

J 

where G pv = SF P ]) + 8Fj$ — 5F$) — 5F$ is the field strength for a non-diagonal subgroup of 
t/(l) 4 . We will continue to ignore this extra field for two reasons. Firstly, it makes the fluctuation 
analysis much more complicated. Secondly, from our condensed matter point of view, we need only 
one gauge field to model electricity and magnetism, and it is not clear to what this other gauge field 
would correspond. In appendix [C] we give more details about this coupling. 
The black hole has metric 

h ds2 = i t- f ^ dt2 + dx ' 2 + dy2 ^ + h%r y ( 34 ) 

and carries both electric and magnetic charge 

F = ha 2 dx A dy + qadz A dt , (35) 
where q, h and a are constants. The function 

f(z) = l + (h 2 + q 2 )z 4 -(l + h 2 + q 2 )z 3 . (36) 
The Hawking temperature, dual magnetic field and charge density of this black hole are [7 
n u 2 2 rp a(3-h 2 -q 2 ) 

B = ha , p= —qa oo , T = . (37) 

47T 
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Note that q 2 + h 2 < 3, whereas p and B are unbounded. We have expressed the charge density in 
terms of the the conductivity of the dual CFT at p = B = 0, 

cr = 2L 2 /kI = V2N 3/2 /6tt . (38) 

In this last equality, we re-expressed the dimensionless L 2 / k 2 in terms of the number N of coincident 
M2 branes [TS]. These relations can be inverted to give a(T, p, B) as the solution of 

3a 4 - AirTa 3 = B 2 + p 2 /a 2 . (39) 

We would now like to consider small fluctuations of the scalar fields (f> and \ about this black 
hole background. It is useful to diagonalise the quadratic Lagrangian (32) first by defining 

k f±^ and ^_ = 4^L. (40) 

The two fields will be dual to two relevant operators 0± . For our dyonic black hole background 

L 4 F^F^ = 2z i (h 2 - q 2 ) , ^e^F^F^ = -Az A hq . (41) 



Substituting into the action (32 1 one obtains the following decoupled actions for the scalar fields 

S± = -^j/ d^xV^g (-\{dM{d^ ± ) + z\q 2 + h 2 )} . (42) 

The equations of motion are therefore 

□V>± = (m 2 ±2z 4 (q 2 + h 2 )/L 2 )i; ± , (43) 

where m 2 L 2 = —2. Note that this mass is above the Breitenlohner-Freedman bound m 2 L 2 > —9/4. 
Depending on the sign of ±, the electromagnetic field provides a potential for the scalar that either 
pushes the scalar away from or towards the horizon. We will assume that the scalar field has the 
following dependence on x, y, and t: <fi ~ e ikx-iuit_ Thus, the equation of motion becomes 

z 4 (j^tf±) - <? z2 i>± + » 2 y V>± " m 2 L 2 iP ± T 2z\q 2 + h 2 )^ ± = , (44) 

where we have defined the dimensionless frequency and wavevector to = tv/a and q = k/a, and also 
/' = df/dz. Near the boundary, z = 0, of this asymptotically AdS± geometry, the scalar has the 
usual behavior 

1>± = z 3 ~ A (A(lu, k) + 0{z 2 )) + z a (B(lu, k) + 0(z 2 )) , (45) 

where A(A - 3) = m 2 L 2 . 

In the standard M2 brane theory, representation theory of 5*0(8) and the superconformal sym- 
metry guarantees that <j) has conformal dimension Aq = 1 while the pseudoscalar \ nas dimension 
Ae> = 2 [T5J H31 H3] . Introducing nonzero B and p, as we have seen, mixes <j> and \- We specu- 
late that in the dyonic black hole background, it is the linear combinations ip± which have definite 
conformal dimension rather than cj> and \ independently. We shall avoid these subtleties by taking 
both 4> and \ to have dimension A<p = 1. The standard M2 brane theory is related to this choice of 
conformal dimension by a deformation by a double trace operator fTJ]. This choice has the added 
benefit of providing a second scalar to analyze with conformal dimension one, and therefore relevant. 

The dimensions Aq = 1 and Aq = 2 both correspond to the mass m 2 L 2 = —2. The fact 
that there are two possible conformal dimensions associated with the same mass has important 
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consequences here. Having chosen the smaller conformal dimension, it is the faster falloff that is 
dual to a source for O in the boundary theory [32] 



Jo = a 3 - Ao A(uj 7 k) , 
whereas the slower falloff gives the expectation value 

(0)= ^a A °(2Ao-3)B(uj,k) 
From these expressions one can read off the Green's function 

(O) = Gg (c,fc)J . 

Thus 



(46) 



(47) 



(48) 



cr B{u,k) 



oo 



(49) 

A/B. Because 



00v -'' 4aA(u,k) 

The standard procedure for scalars of mass m 2 L 2 = —2 would have yielded Gq 
our operators have the smaller conformal dimension, we get the result (49). 

Our job is to determine the dimcnsionlcss function F appearing in the scaling relation (13 1. We 
will need to compute two scaling functions F±, one for each scalar field ip±. An observation that 
simplifies our task is that the differential equation ( 44 ) depends on B and p only in the combination 
h 2 + q 2 , or equivalently B 2 + p 2 /a 2 . Thus, F is a function of B 2 + p 2 /a 2 onljj^] Using the fact that 
for the dyonic black hole [Zl [9] 



0-ocr 



(1 



P 



e 

2 ■ 



and that for the M2 theory (12 1 becomes 



1 



= — F 
T 



1 

J^2 



l B 2 + p 2 /a 2 

the expressions in the preceding paragraphs imply that 



(50) 



(51) 



F 



sT 



167r 2 (e + P) xo^o 



lim / dq q Im 



wA(m, q) 



(52) 



It is interesting to recall from |9] that the dc conductivity for the M2 brane, with nonvanishing B 
and p is given by 

_ isT) 2 



allowing us to rewrite ( 52 1 as 
1 



F 



rp2 



B 2 + P 2 lo 2 = 



16tt 2 



<70 ro->0 



lim / dq q Im 



B(to,q) 
ro^l(ro, q) 



(53) 



(54) 



The physical meaning of this formula is not clear, because in general there's no particular reason 
for the impurity relaxation timescale to be related to the electrical conductivity of the CFT in the 
absence of impurities. It remains for us to solve for A(ro,q) and B(to,q) numerically. 

3 Amusingly, this dependence only on the combination h 2 + q 2 remains true in the full M2 brane theory, in which 
one keeps the second gauge field to which cj> and x couple. See appendix [C| 
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4.1 Comment on discrete symmetries 

Before turning to computations, we should comment that we have found two relevant operators 0±, 
both with dimension Aq ± = 1. Both operators preserve the global C/(l) symmetry and in principle, 
they could both be sourced by impurities. However, these operators transform differently under the 
discrete symmetries charge conjugation C, parity V, and time reversal T. 

Usually we make a sharp distinction between the way operators transform and the way states 
transform. The situation here is complicated by the fact that our fields ip± dual to the 0±, and 



hence the 0± themselves, are defined in (40 1 in terms of the state and in particular the choice of 
background magnetic field B and charge density p. In studying ip±, when we act with C, V, and T, 
we will act on both the operator and the underlying state. 

From the gauged Af = 8 supergravity construction we know that <fi and x are both real fields 
and thus do not transform under C. At the same time is a scalar while x is a pseudoscalar. Thus 
we have V(4>) = <j) and V(x) = ~X- From the CPT theorem, we must have under time reversal that 
T{4>) = 4> while T( X ) = -X- 

The transformation properties of the electromagnetic field under C, V, and T are well known. In 
our case, the strength of the electric field of the black hole, which is dual to the charge density p on 
the boundary, is characterized by the value of q. The magnetic field strength B is fixed by the value 
of h. Under parity, we know the electric field is a vector while the magnetic field is a pseudovector. 
Thus we have V(q) = —q while V(h) = h. C flips the sign of both q and h. Finally, the magnetic 
field famously breaks time reversal symmetry, and so T(h) — —h while T(q) = q. 

Assembling these facts about q, h, <fi and x together we find that 

Cty+) = -</>+ , T(iJj+) = i>+ , T(i,+ ) = -V+ , (55) 

The subscript of rp± thus corresponds to the eigenvalue of the field under V . In adding impurities, 
we are faced with a choice: we can let the impurities break T but not V or V but not T. In the 
next section, we will see the behavior of the scattering time T lmp for the 0+ impurities is markedly 
different from the 0_ impurities. 

Note that 0+ transforms the same way under the discrete symmetries as the magnetic field while 
O- transforms as the electric field. It is thus tempting to call the 0+ impurities magnetic and the 
0_ impurities electric. The underlying scalar fields 4> and x are °f course neutral. The nontrivial 
transformation of ip± under C comes from the dressing by the dyonic black hole background, or 
equivalcntly, the presence of magnetic field and charge density in the field theory. 

5 Numerical results for l/ri mp 

The prescription for computing the retarded Green's function |28j is that we solve the equation for 



our scalar field in the black hole background ( 44 1 with ingoing boundary conditions at the black 
hole horizon z = 1, 

i> ~ (l - 2 )»">/c>W-3) _ (56) 

It is convenient to impose these boundary conditions by defining a new function S(z) such that 

i()(z) = exp (itv J J^ dx ^j S ( z ) ■ ( 57 ) 
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Thus S(z) is nicely behaved at the horizon and we can impose S(l) — 1. It is straightforward to 
integrate the equation for S(z) numerically from the horizon to the boundary at z = 0. Near the 
boundary one can read off the value of the ratio B/A that we need to compute the scaling function 
F. 



Given B/A, we then need to compute the integral (52 1 over k. This integral may not converge, 
and indeed we discuss below a large p and B regime where it does not for the ^_ scalar field. 
However, the divergence comes from bad behavior at finite k. At large k we expect the imaginary 
part of the Green's function to be exponentially damped. Appendix [D] presents a WKB analysis 
that demonstrates this damping. Qualitatively, we may understand the behavior from the bulk as 
follows. At large k, the Green's function is dominated by a space-like geodesic which stays close 
to the boundary of the space-time and far from the black hole horizon; it is the horizon that is 
responsible for producing dissipative effects and a nonzero contribution to the imaginary part of the 
Green's function. 

We begin by considering the two point function corresponding to the scalar field ip+, that is, 
magnetic impurities. The result for F + is plotted in figure IX^u. As the combination B 2 + p 2 /u 2 
is increased, the inverse scattering time becomes smaller and smaller. There is also a competing 
physical process that we described in section 3.2 above: the cyclotron resonance. If we keep B/T 2 



small (h <C 1), then we can trust our hydrodynamic approximation of the location of the cyclotron 
pole at lu c — vy, given in (30 1. In figure |2j we plot the dimensionless j/T against p/T 2 <jQ for an 



appropriately small B/T 2 = 1. We see that at large p, the damping due to the cyclotron resonance 
also gets small while the frequency changes more slowly. The reduction in damping both from 
impurities and the B field at large p may make the hydrodynamic cyclotron mode easier to observe 
experimentally than the results in [5] suggested, as there the dependence of T; mp on B and p was 
not considered. 

We also would like to consider the scaling function F- for the ip- scalar fields. The result is 
plotted in figure [T]d, and is somewhat dramatic. We define the related dimensionless quantities 

Q^^WT? and QeI^+^K 2 • (58) 

We find that F_ increases monotonically as a function of Q or Q and at a finite value = 0.76654, 
which corresponds to = 20.80, l/r imp diverges. This divergence reveals some interesting physics 
in the 'pure' CFT which we will return to shortly. Note that we cannot trust our memory function 
formalism if l/r imp gets too large. Higher order and nonperturbative effects are important in this 
regime. 

From a gravitational standpoint, we may attempt a qualitative understanding of the decrease 



(increase) of l/ri mp with p and B for the ip+ (tp—) field. From (44 1, the scalar ip± experiences a 
potential of the form ±2z 4 (/i 2 + q 2 ). Thus for the potential repels the scalar from the black hole 
horizon while for the potential has the opposite effect. The black hole horizon is responsible 
for dissipation in the system, and we expect the more the probability distribution of the scalar is 
peaked near the horizon, the larger the imaginary part of the scalar Green's function will be and 
hence the larger 1/Ti mp will be. The decrease in F± with B and p for and the corresponding 
increase for tp_ in figures [l^, and [TJd are thus qualitatively understood. 

5.1 A black hole instability 

The divergence in F— is caused by an underlying instability in the truncated M2 brane theory at 
large Q. Numerically we have found that the divergence in F— appears to be caused by a diffusion 
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Figure 2: a) The real part of the cyclotron frequency and b) momentum relaxation due to the 
cyclotron resonance at B = T 2 . 

type pole in the scalar Green's function moving onto the real k axis. Near the critical value Q*, the 
pole has the approximate form 

G °° ~ C(M)-q2 + 6itt)/5+... ' (59) 
The ellipsis denotes higher order terms in u> and k and C(h, q) — 1.3911(Q — Q*). In figures]^ and 
b we have plotted the location of the pole as a function of Q. Note that in both plots, for k 2 and 
Imw, the curves reach the real axis at = 20.80. At this point, two changes occur simultaneously. 
First, l/r; m p diverges because, at u> = 0, there is a (double) pole on the real k axis that is integrated 



over in (52 1. Second, at k = 0, the imaginary part of ui switches sign leading to exponential growth 
in the scalar mode at long times. 

This exponential growth in time is indicative of an underlying instability in the CFT when 
Q > Q*. Indeed, the retarded Green's function cannot have poles in the upper half lo plane, 
indicating that for Q > Q* we are no longer in the correct vacuum. To evaluate the impurity 
time for Q > Q*, we would need to find a different stable supergravity background in which the 

field was at a local minimum. This instability illustrates the dangers in this particular black 
hole background of working in a consistent truncation in which we set -0_ = 0; had we ignored 
fluctuations in we would have missed this instability. 
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Figure 3: The location of the smallest pole in Gq_ _ as a function of B and p: a) The pole in the 
complex k plane for uj = 0. Note k 2 is real, b) The pole in the complex u) plane for k — 0. Note 
Re(w) = 0. 



This instability in the theory is interesting from a purely gravitational perspective. First note 
that poles in the Green's function (49 1 occur when the denominator A(lj, k) vanishes. This vanishing 
means there exist on shell modes with the behaviour ip± = z + C(z 3 ) near the conformal boundary 
z = 0. Due to the nonstandard boundary conditions for the tp± fields, these modes are physical. For 
the critical value Q — Q*, the physical mode is static, with u) = 0. This 'threshold' mode separates 
oscillating modes, with Q < and Imu < 0, from exponentially growing modes, with Q > and 
Imw > 0. For Q > there is an on shell classical instability of the black hole. 

A classical instability of the full finite temperature M2 brane was discovered in [T51 H7J- To 
see that instability, it is necessary to retain more than one Maxwell field in the bulk, whereas we 
have truncated to a diagonal combination. Converting to our variables, the Gubser-Mitra instability 
occurs with h = and at qcm = 1, corresponding to the value pgm / pqT 2 k, 39.5. This charge 
density is larger than necessary for the instability we have described. Given that the two instabilities 
occur in different theories, one might be cautious about attaching much meaning to the relative values 
of the critical charge^] 

The classical Gubser-Mitra instability occurs at precisely the same value of the charge density 
at which the charged black hole with four U(l) gauge fields becomes thermodynamically unstable. 
This coincidence led Gubser and Mitra to the 'correlated stability' conjecture, that for translationally 
invariant horizons, thermodynamic and classical dynamic instabilities should always coincide. 

In contrast, our black hole with only one U(l) gauge field is always thermodynamically stable 
[TJ; it only becomes thermodynamically unstable when embedded into the full M2 brane theory, 
with four independent charges. Given this dynamical instability without a corresponding thermo- 
dynamic instability, we have a counterexample to the correlated stability conjecture. The physics 
underlying our counterexample appears to be very similar to the counterexamples discussed in [T5] . 
The mismatch between thermodynamic and classical instability is possible because the instability is 
due to scalar fields that are not associated to conserved charges. The instability indicates a phase 
transition to an ordered phase, characterised by a condensate for the operator C_ dual to the field 



4 In fact, the analysis in |16l I17| also employs an inconsistent truncation. Although that work does include the 
gauge field G, it does not include the pseudoscalar \, which is sourced by G. See appendix[C] For an investigation of 
threshold modes u) = k = 0, this xGxy coupling will vanish and so should not affect the location of the Gubser-Mitra 
instability. Note also, 16, T7] work with conformal dimension A = 2 for the operator dual to 0, so their dual theory 
is also related to the standard M2 brane theory by a double trace deformation. 
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%))—. As noted in [13], it seems appropriate to widen the notion of thermodynamic instability of black 
holes in these cases to allow for variations of the value of the scalar field at the conformal boundary. 
Because this phase transition is associated with breaking an underlying Z2 symmetry, ip_ — > —ip—, 
of the truncated supergravity action, the transition is likely to be second order. 



6 The Nernst effect in the truncated M2 brane theory 

Place a system in a thermal gradient and a perpendicular magnetic field. The Nernst effect is the 
observation of an electric field that is created orthogonal to both the thermal gradient and the 
background magnetic field. The Nernst coefficient is the electric field generated per unit of thermal 
gradient and magnetic field 

N = ^— — . (60) 
BWT y ' 

Recent experimental interest in the Nernst effect in superconductors was sparked by the obser- 
vation by Ong et al. of an anomalously large Nernst effect in the pseudogap region of the high-T c 
superconductors |33j . This observation was interpreted as signalling the presence of vortices, and 
hence the persistence of a phase disordered Cooper pair condensate, above the transition tempera- 
ture. The Nernst effect has since been used to probe the presence of long lived Cooper pair fluctu- 
ations in conventional superconductors |34j and also to capture the proximity of a Mott insulating 
phase in certain organic superconductors [19] . 

The Nernst effect depends strongly on the impurity relaxation time. Therefore, the B and p 
dependence of T; mp that we explored above may have measurable consequences for observations of 
the Nernst effect in various types of superconductors. In this section, we explore the dependence of 
the Nernst effect on rj mp . We argue that there are qualitative similarities between aspects of the 
T and B dependence of the Nernst effect observed in organic superconductors [19] and the Nernst 
effect we predict at p = 0. In the underdoped cuprate superconductors [33J [35], f° r a range of 
typical values of the doping, we find that the relatively large size of p compared with experimentally 
accessible values of B make Tj mp relatively insensitive to B. The exception to this last statement is 
very close to the insulating state at doping xj — 1/8, where p = is again appropriate. 

A hydrodynamic approach to the Nernst effect was developed in ref. [3J. That work led to the 
following formula for the Nernst coefficient in a relativistic theory with 'speed of light' v. For the 
next couple of formulae only we shall not set v — 1, as we have been doing in this paper, and shall 
also show explicitly the unit electric charge e. It was found that 

v 2 ( l/n 



T V(^/7+l/r im p) 2 +c c 2 
Restoring the v and e dependence of the cyclotron pole, we have 



r • (6i) 



eBpv 2 aB 2 v 2 
uj c = — — - and 7 = — — . (62) 



In ref. [5] the hydrodynamic expression (61 1 was used to reproduce some features of the experi- 
mentally observed Nernst effect in the cuprates. However, in the absence of a microscopic theory, the 
simplifying assumption was made that the functions of state, the conductivity a and the impurity 
scattering time Ti mp did not depend significantly on the dimensionless ratios p/T 2 and B/T 2 . For 
the M2 brane theory, <r(p/T 2 , B/T 2 ) was computed in ref. [5J, and in this paper we have computed 
Timp- Thus we have all the ingredients necessary for an exact computation of the Nernst coefficient 
at a quantum critical point. 
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The Nernst effect is a direct probe of the impurity relaxation time. We see in (61 1 that a nonzero 
Ti mp is necessary for the Nernst coefficient to be finite. This connection becomes particularly striking 
in the limit of vanishing charge density, p = 0, where (61 1 becomes simply 

2 2 

at v T imp V (f\X\ 

T ~ V 2 F(B/T 2 ) ' ( ' 



where in the second step we used (51 1 for the M2 brane theory. One nice feature of this formula 
is that the dependence on 7 has dropped out, and with it the uncertainties of order V 2 B 2 that we 
discussed at the end of section 13.21 



Using our numerical results from the previous section, it is straightforward to plot the Nernst 
coefficient for the truncated M2 brane theory as a function of B and T, which is how the experimental 
data is often presented (TSJ [33 13S] • In figure]^ we plot the Nernst coefficient in the case of vanishing 
charge density (63). We have taken the impurity potential to couple to the 0+ operator. If we had 
taken the C_ operator, the dark and light regions of the plot would have been interchanged. For 
ease of comparison with experiments, we should consider plotting in SI units. In order to convert 
B/T 2 from our 'natural units', with v — ks = h = e = 1, one must first introduce the velocity 
v = <!> v m/s, and then note that 

= 8.8 x HT^i ■ ( 64 ) 



(kelvin) 2 ' v ehv 2 

We now need to specify the number A convenient choice is = 3400, as this makes the 
conversion factor 8.8 x 10~ 8 Q 2 = 1. Curiously, this appears to be in a physically reasonable ballpark, 
at least for the cuprate superconductors. This value corresponds to hv — 22meVA, whereas the 
value estimated in [5] for underdoped LSCO was 47meVA and the characteristic velocity observed 
in YBCO was 35meVA [37]. While all plots in this section are in natural units, in figure^, for 
concreteness only, we have assumed hv = 22meVA and put kelvin and tesla on the axis labels. 

In figure [4^, the range of temperatures and magnetic fields has been constrained by the validity 
of the hydrodynamic approximation, which for the M2 brane theory requires B/T 2 < 10 [9 . Outside 
this regime one does not expect (61 ) to hold. Although the M2 brane theory can be studied beyond 
the hydrodynamic regime [9], in that case our arguments in section 3.2 for combining the effect of 
impurities and the magnetic field do not hold. 

The qualitative form of figure^ is determined from the fact that N ~ 1/ F(B /T 2 ), with F a 
monotonically decreasing function. This scaling will be a universal feature of the Nernst effect at 
p = in the vicinity of a quantum critical point. The scaling near a general quantum critical point 
will be 

N ~F<B]Wy ^ 

We illustrate this general form in figure[4j3, taking F(x) = 1/ (1+x 2 ) for concreteness. It is interesting 
to note that this qualitative form appears to arise in Fig. 3b of ref. [TH] at low magnetic fields and 
just above the superconducting phase transition temperature, T c . The authors of ref. [19 study 
organic superconductors, and their Fig. 3b is, like our figure [4] a plot of N(T, B). Their system also 
has zero charge density, as it is not doped. The region of interest (low B, just above T c ) is precisely 
that for which 19J propose that strongly correlated electron physics is important, as are the effects 
of a nearby Mott transition. 

We can try to compute the Nernst signal with a finite charge density using the formula (61 1. As a 
benchmark, one can compare with experimental results for the underdoped cuprate superconductors. 
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Figure 4: (a) The Nernst signal N for the truncated M2 brane theory as a function of B and T, 
with vanishing charge density p = 0. The impurity potential is coupled to O+. Lighter denotes a 
larger Nernst coefficient, (b) The Nernst signal in the vicinity of a general quantum critical point 
with p = and Aq = 1. The shading is logarithmically spaced. 

It was argued in [8] that the charge density was given by the difference in the doping from the 
commensurate insulating state at xj = 1 /8, divided by the area of a unit lattice cell. For underdoped 
LSCO, at x — xi — —0.025 for example, the lattice constant is a = 3.78 A, leading to 

p= e ^ x ~ Xl ^ ^ 0.028-^. (66) 

The conductivity at the critical point can be estimated to be [5H] 

4 e 2 Q2 

ao w — w 1.6 x 10- 4 — , (67) 
h J s 



175tesla. (68) 



allowing us to obtain 

P_ 

We have expressed the result in tesla, so that it may be compared to the value of the background 
magnetic field. The value of 175 tesla is significantly larger than the typical magnetic fields applied 
in experiments measuring the Nernst effect in the cuprate superconductors. In our expression for 
the Nernst coefficient, the magnetic field B only appears in the combination B 2 + p 2 /<7q. Therefore 



the large value of the charge density in (68 1 swamps out the B dependence, and the resulting plot 



of N{T,B) shows simply vertical lines. We might note however that both (66) and (67 1 are only 



tentative identifications, so the result (68) could change substantially. In particular, the conductivity 
( 67 ) assumes the putative nearby superconductor-insulator quantum phase transition in the cuprates 
is in the same universality class as the films reviewed in [38 . 

It was suggested in [8j that the quantity that should be compared with measurements of the 
Nernst effect in the cuprates is not N, but rather the off diagonal thermoelectric coefficient a xy . This 
was to isolate the contribution of critical superconducting fluctuations from non-critical fermionic 
contributions to the condutivity. This coefficient gives the electrical current generated by an applied 
thermal gradient: J x = —a xy d y T. It is related to the Nernst response by the electrical conductivity, 
roughly a xy ~ a xx BN. Here a xx is the full conductivity of the system, not just that of the critical 
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fluctuations. The thermoelectric coefficient has a B dependence that is not swamped out by a large 
charge density [8] , and indeed reproduces various features of the experimental data for N(T, B) in 
the cuprates. The large charge density implies that if we use our M2 brane expressions for the 
impurity timescale Ti mp (p, B, T), we will not introduce any extra dependence on B into the previous 
results of [5] • In any case, the high critical temperature of the cuprates combined with the value for 
the velocity v, discussed above, implies that B /T 2 , p/T 2 <C 1 over the region of interest. Therefore 
it is consistent to neglect the B and p dependence of Ti mp and the e + P for these systems. 

In figure [5] we have plotted the Nernst coefficient for the truncated M2 brane theory with p/<Jo = 
8.3 x 1CP 5 and p/ao = 10, in each case with V = 0.1. We have again used 0+ as the operator 
coupling to the impurity potential. The plots are significantly different and less universal than the 
vanishing charge density case, of figure [4] As p becomes nonzero, a local maximum appears at small 
temperatures and moves rapidly to the right as the charge density is increased. Note that unlike 
the p = case we considered previously, plot [5^, in particular is vulnerable to the unknown V 2 B 2 
corrections to 7 that we mentioned in section |3.2| above. Its qualitative form however should be 
correct, as it needs to interpolate between figures [3] and [5b. 
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Figure 5: The Nernst signal N for the truncated M2 brane theory as a function of B and T, with 
(a) charge density p/ao — 8.3 x 10~ 5 and (b) charge density p/a n = 10, and V = 0.1 in both cases. 
The impurity potential is coupled to C+. Lighter denotes a larger Nernst coefficient. 

What we learn from these plots is that (with V = 0.1) a very small p/<Jo ~ 10 -5 is already 
sufficient to produce significant deviations from the p = result. This suggests that one would have 
to tune very close to the insulating doping of xj = 1 /8 in order to see the scaling p = behaviour 



of equation (651. However, it is curious that the observations in [35 do appear to show similar 



behaviour at low B and with T just above T c 



7 Discussion 



7.1 Summary of results for the impure CFT 



In this paper we have derived a general formula ( 12 1 for the impurity relaxation timescale at a 
quantum critical point. Some of our steps in section 3 were adapted from earlier works [25, 26, 27 . 
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We then used our formula to obtain T lmp for the truncated M2 brane theory. We found that there 
were two relevant operators 0± that were neutral under a global U(l) symmetry. This led to the 
following results 

• For the 'magnetic' operator dual to the scalar field ip+, increasing the magnetic field B and 
charge density p suppresses momentum relaxation due to impurities, see figure |TJl. At small 
B and sufficiently large p, this implies that the cyclotron resonance described in [5J [S] may 
be easier to observe experimentally than previously estimated because the impurity damping 
effects become smaller. 

• For the operator dual to the scalar field in contrast, the inverse relaxation time increases 
when the charge density p and magnetic field B are increased. Beyond the critical value 
Q = \J B 2 + p 2 /<Jq ~ 2 IT 2 , the inverse relaxation time actually diverges. See figure^. We 
saw that the divergence is associated with a new instability in the underlying finite temperature 
CFT. The scalar field ip- develops an unstable mode when the Q gets too large. This instability 
is dual to a phase transition to an ordered state in which 0_ condenses. 

• The Nernst coefficient for the truncated M2 brane theory with vanishing charge density, p = 0, 
is shown in figure [4] with impurity relaxation due to the 0+ operator. We noted that our 
dependence on T and B is qualitatively similar to the results for low magnetic fields and for 
temperature close to the critical temperature observed in various types of superconductors 

[35J EM] • This comparison makes most sense for the data in [19 , in which the system is not 
doped. It is interesting to note, however, the similarity in the Nernst response between the 
different superconductors in this regime. 



7.2 Some directions for future work 

Our results are perturbative in the strength of the impurity potential V and valid at high tem- 
peratures. We furthermore used the memory function formalism, which requires l/ri mp not to be 
too large. Various techniques have been developed in condensed matter theory for treating random 
disorder more exactly. It would be very interesting to see if these can be applied to the M2 brane 
theory. 

Our treatment of the M2 brane theory was not completely satisfactory. We made an inconsistent 
truncation in order to avoid dealing with more than one U(l) gauge field. We also made a noncon- 
vcntional choice for the conformal dimension of our supergravity mode x- One clear project for the 
future is to repeat the calculations presented here including the extra U(l) gauge field and using 
the standard choice of conformal dimension for \. 

But it would also be interesting to take a more general approach. While in our M2 brane 
theory, the operators we considered had conformal dimension Aq = 1, in general the spectrum of 



operator dimensions will depend on the CFT. Our bulk equation of motion (44 1 for the scalar field 
4> dual to the operator O can be generalised to allow the scalar field to have an arbitrary mass. 
This is dual to allowing Aq to vary. One could imagine using the AdS/CFT correspondence more 
'phenomenologically' to see how the effects of impurity relaxation depend on the dimension of O. 

One result of this paper is the instability of the truncated M2 brane theory at \J B 2 + p 2 /cr^ ~ 
21T 2 . We noted that this instability is a counterexample to the original formulation of the correlated 
stability conjecture, as the black hole is thcrmodynamically stable at this value of the charge density 
and magnetic field. As with the previously studied counterexamples |18j . the mismatch between 
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dynamics and thermodynamics occurs because there is no conserved charge associated with the 
scalar fields. The instability presumably indicates the existence of new dyonic black hole solutions 
in which both the scalar (j> and pseudoscalar \ are nonvanishing. It would be interesting to find 
these new solutions and understand their implications for the dual theory. 

In combining the effects of impurities and an external magnetic field, we were constrained to 
work in a regime where B/T 2 was not too large. It would be nice to relax this requirement. An 
approach via Ward identities looks promising. 

We have noted that the Nernst effect in systems with vanishing charge density can be a very 
clean indicator of the proximity of a quantum critical point, via the scaling relation (65). Perhaps 
this scaling can be searched for in future measurements of the Nernst effect. 
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A Relating Gp V to a conductivity 



The connection between G^p-p(u>,0) and 9 is established as follows. The momentum conductivity is 
defined by 

(7>}=-0(w,O)VT. (69) 

The momentum density V can be sourced by fluctuating the background spacetime metric. Specifi- 
cally, by definition there is a coupling between metric fluctuations and the stress tensor in the action, 
SS= §/ cPxT^Sg^. (We are treating T^ v as a tensor density rather than a tensor field and hence 
have absorbed a factor of y/—g.) Without loss of generality by rotation invariance, we can restrict 
to a fluctuation only in the spatial direction Sgoi- Choosing n 3 = 5ij, linear response theory then 
implies 

(V) = G^ v (io,0)5g 0j n 3 . (70) 

The remaining step is to show that Sgoi is gauge equivalent to a thermal gradient. Recall that under a 
diffcomorphism generated by the vector £ a , metric perturbations transform as Sg a i, = d a £b-\-db£ a . We 
choose £j = and di£o — —Sgoi. After this gauge transformation Sgoi vanishes and diSgoo — 2iu>5goi, 
assuming the fluctuations have a time dependence of the form e~ luJt . Recall that the Euclidean time 
direction is periodic with period l/T. It is convenient to fix the period and take goo = 1/T 2 . A shift 
in temperature thus implies dtgoo = —2diT/T 3 — — 2goodiT/T, to leading order in d{T. Putting 
these formulae together leads to the intermediate result 

S90i = -^r- (71) 

We have to be careful in interpreting this result because we rescaled time at an intermediate step. 
Moreover, we rescaled goo without rescaling uj. We now need to rescale time back to the lab frame, 
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where time has period l/T in the Euclidean direction, being careful not to rcscale lu: 

djT 
iojT 



Sgoi 



Thus we find from (69 1 and (70 1 that 



(72) 



(73) 



as claimed in the main text. As we noted, &p V (oj,0) needs to be corrected by a contact term to 
give a sensible oj — > limit for the conductivity. 



B Hydrodynamics and the memory method 

In this appendix we show that the results obtained from hydrodynamics in [9] and [8] are completely 
consistent with the memory function methods. Consider Gjj = oja + for the choice n,i = (1, — i)/y2. 
We will henceforth suppress the superscript R since it is clear we work only with retarded Green's 
functions in this appendix. In [9 , for a translationally invariant theory in the hydrodynamic limit, 
it was argued that this holomorphic conductivity takes the form 

. w + iu*/'y + u c 

a + = ig . (74) 

oj + v~f — tu c 



Using ( 25 1 , we find that 



G vv {oj) = {e + P) 7. Uc . (75) 

OJ + Vf — 0J C 

As expected, Gp>v(fi) = e + P = \o- Moreover, this result is entirely consistent with the memory 
function formalism, 

LuG P p{uj) . 

hm = vy — uj c . (76) 

b^o xo 

This limit is B/T 2 — > 0. The frequency oj/T is kept small but fixed. 

In [8] , it was shown that an impurity scattering time could be introduced by making the replace- 
ment oj — > lu + i/r in the conductivities. This followed from hydrodynamics plus a relaxation of 
momentum imposed by hand. This substitution yields 

, . u + i/T + iu 2 /j + u c 
cr, = ia —. ■ . (77) 

T OJ + l/T + IJ — LO c 

which agrees with the combination a xy + ia xx from 0. We use the superscript ' to indicate a Green's 
function or conductivity in the presence of impurities. The prescription is trickier when applied to 
G-p-p. The relevant transport coefficient we will call 9 + is related via ujT9 + = G-p-p(uj) — xo- Thus 

6'T = (78) 

OJ + l/T + IJ — Ld c 



G vp {oj) = xo . , , , ,_ , ,., — — • (79) 



which then implies that 

i/r + if — oj c 
' uj + i/r + ij — lu c 

Again, the memory function method works with this modified G' vv yielding a pole at the appropriate 
lu = —i/r — if + oj c . 
Let's now compute 

„9/ / / \ ,„,, (oj + i/T)(ij — uj c ) . n . 

B 2 (a' + (u) - (7+ = Xo , 1 , .; '2 . Cj s , (80) 

{OJ + l/T + 17 — LU C ) 
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and check the extent to which (251 continues to hold in the small 1/r and B limit: 



lim — (g' vv {u)-~{G , J j{u>)-Gjj{0))) =i 1 - U c + i/T-{i 1 -U c )=i/T . (,sl) 
B A— XO V 



Indeed, that's exactly the result we were looking for, suggesting (29) is correct. We have made only 
an indirect assumption about the B dependence of 1/r. Our hydrodynamic framework is valid only 
up to 0(B 2 ). Thus while (22 1 appears to capture the full B dependence of 1/r, in this hydrodynamic 



context, we can trust the answer only up to 0(B 2 ). 

C Action for fluctuations without truncation 

In this appendix we give, for possible future use, the untruncated action for perturbations of the 



dyonic black hole background ( 34 1 and ( 35 1 



From appendix F of [30] one can check that the full action for 3 scalars, 3 pseudoscalars and 4 
gauge fields may be consistently truncated for linearised perturbations of the dyonic black hole to 
one scalar, <fi, one pseudoscalar, x, and one gauge field G: 



S = - 

2k 



K f dtx^g (R - \ {(d^)(d^) + (d^xWx)} + 2£- 2 [3 + \(<t> 2 + x 2 )] 
-L 2 [1 + |(0 2 - X 2 )] F^F^ + \L 2 ^ lvpa F^ u F pa 
-tfG^GT - 2L 2 <j>F^G^ + L 2 X e^ xp F^G Xp ) . (82) 



We have used the same notation as in the main text. In particular G pv — 5Fj£) +5F$ —5F^) —5Fj$ , 
whereas the background gauge field is F pv = F$ + F$ + F^J + F$ . 



Once again it is useful to introduce the fields ip± defined in ( 40 ) . Evaluated on the dyonic black 
hole background, the action becomes 

S =-^2/ rf W=S (-|[(^+)(^+) + (^-)(^_)]+i- 2 (^ + V-) (83) 

-L- 2 z\h 2 + q 2 )(i; 2 + -i> 2 _)~ L 2 G pv G^ (84) 



-4L- 2 ^ 

i \ -■ 



V/i 2 + q 2 (a^G tz + ip+Gxy) ) . (85) 



In this expression we note that once again the background magnetic and electric charges only appear 
in the combination h 2 + q 2 . It should be possible to decouple the equations of motion following from 
this action and find the momentum relaxation in the full theory. Solving the full equations would 
also be of interest in terms of revisiting the computations in |16l 117] . 

D Exponential falloff at large k 

The exponential falloff of the imaginary part of the Green's function with large spatial momentum 



k can be seen from a WKB analysis. To perform a WKB analysis of the differential equation (44 1, 
we first transform it into Schrodinger form. (A very similar analysis was performed in an appendix 
of ref. [28 .) To that end, we define a new wavefunction <i>± such that &±(p) = tf)±(z)/z and a new 
radial coordinate p(z) such that d z p = 1//. With these new definitions, the differential equation 
becomes 

(-d 2 p + V(p)-tv 2 )<i> ± (p) = (86) 
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where 

V(p) = ^(m 2 L 2 + qV + 2/ - zdj ± 2z\q 2 + h 2 )) . (87) 
In terms of the new radial variable, we can take 




Thus, the horizon at Zh = 1 has been moved to ph — > oo. The boundary at Zb = remains at p^ = 0. 

Note that for m 2 L 2 = —2 we have V(0) = q 2 . We are interested in the case of spacelike 
momentum, k 3> to, for which there are no classical turning points. The wavefunction is exponentially 
damped everywhere. The imaginary part of the retarded Green's function will be proportional to 
the probability for the particle to reach the horizon at Zh = 1: 




12 3 4 

k/a 



Figure 6: The points are numerically determined values of the log of the imaginary part of 
G@ o (0,k) as a function of k for q 2 + h 2 = 1/4 and the ip+ scalar field. More specifically, we 
plot ln(lim ro ^o B/wA). The line is a best fit 1.15 — 3.15q. For these values of q and h, WKB gives 
us a = 2.86. 

Given that k ^> oj, we can approximately evaluate this integral to give ImGg ~ exp(— ah/ a) 
where we have defined the coefficient 

dz ,„„. 

7m ■ m 

In other words, the WKB analysis is telling us that Im Gq is exponentially damped for large enough 
k/a. In the case where h = q = 0, we may evaluate the integral analytically to find 

dz -^W 3 ^ 1.402. (91) 



\/r^3 r(5/6) 



In order for our integral (52 1 over q = k/a of the imaginary part of the Green's function to converge, 
this exponential damping is important. We observed it numerically for more general values of h and 
q. In Figure [6j we plot the log of the imaginary part of Gq O (0, k) as a function of k/a and indeed 
see linear behavior at large k. 
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